REPORT  DOCUMENTATION  PAGE 

form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  colieaion  of  information.'  Send  comments  regarding  this  burden  estimate  or  any  other  aspert  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Proiea  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

31  Mar  95 

3.  REPORT  TYPE  AND  DATES  COVERED  | 

final  30  Sep  93  through  31  Mar  9p 

4.  TITLE  AND  SUBTITLE 

Fracture  Mechanics  of  Advanced  Composites 

(DEPSCoR  92) 

5.  FUNDING  NUMBERS 

G  F49620-93-1-0550 

6.  AUTHOR(S) 

Michael  C.  Larson 

AFOSR-TR-95 

ObVA 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

1  Department  of  Mechanical  Engineering 
;  Tulane  University 
;  New  Orleans,  LA  70118-5674 

r 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

[  AFOSR/NA 

i  110  Duncan  Avenue  Suite  B115 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

1 

Bolling  AFB,  DC  20332-0001 


11.  SUPPLEMENTARY  NOTES 


19951027  098 


r 

.  - .  ...  - 

_ nfrr: _ 

: 

i  12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

I  -  . 

(^ELECTEI 

NOV  0  3  1995 ; 

^  ..  I 

L_ 

)ISTRIBUTlON  CODE 

DTIG  QUALITY  INSPECTED  5 
This  report  documents  instrumentation  obtained  under  a  DEPSCoR  grant  to  support  research  in 
the  mechanics  of  advanced  composite  materials.  While  the  equipment  will  support  a  broad  range  of 
fresearch  topics,  it  has  already  spawned  studies  into  micromechanical  fracture  of  advanced  composite 
toateri^s  and  nondestructive  evaluation  of  cracked  bodies.  Specifically,  a  laser  interferometry  system 
(including  optical  and  vibration  isolation  components  with  image  processing  hardware  and  software),  and 
|a  tensile  substage  for  an  electron  microscope,  form  the  core  of  an  experimental  program  which 
^complements  analytical  and  computational  investigations. 

:  Progress  is  being  made  in  two  particular  areas.  The  first  of  these  is  the  development  of  a  coupled 

iexperimental  and  computatipnal  methodology  for  solving  inverse  fracture  problems  using  laser 
^olographic  interferometry.  The  methodology  can  be  utilized  in  the  fracture  mechanics  of  homogeneous 
sand  heterogeneous  media  to  determine  the  location,  size,  and  geometry  of  an  embedded  crack  given  the 
imeasured  displacement  history  of  a  free  surface.  The  technique  will  be  valuable  as  a  research  tool  for 
-those  interested  in  tracking  the  development  and  propagation  of  cracks  inside  test  specimens.  The  second 
program  involves  examining  microscale  sliding  and  debonding  at  fiber-matrix  interfaces.  This  study 


1  14.  SUBJECT  TERMS 

1  Fracture,  Inverse  Problem,  Holographic  Interferometry 

1 

a _ 

15.  NUMBER  OF  PAGES 

13 

16.  PRICE  CODE  1 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL  i 

NSN  75^0-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std  Z39-18 


HNAL  REPORT 


for  the  period  September  30, 1993  through  March  31, 1995 
Grant  Number  F49620-93- 1-0550 


Fracture  Mechanics  of  Advanced  Composites 
(DEPSCOR  92) 


t  ;  ( // 

Michael  C>Larson 
Department  of  Mechanical  Engineering 
Tulane  University 
New  Orleans,  LA  70118 
voice:  (504)865-5775 
fax:(504)865-5345 
email:  larson@me.tulane.edu 


submitted  to: 
Contracting  Officer 
AFOSR/PKA 
1 10  Duncan  Avenue 
Suite  115 

Bolling  AFB,  DC  20332-0001 


Contents 


I.  Overview . 3 

n.  Equipment  Acquired . 4 

in.  Appendix . 6 

Publications  describing  investigations  which  utilize  the  equipment 
obtained  under  this  grant. 


2 


Overview 


1.  Overview 


The  equipment  obtained  under  this  grant  is  supporting  broad  research  in  applied 
mechanics  and  has  already  spawned  studies  into  micromechanical  fracture  of  advanced  composite 
materials  and  nondestructive  evaluation  of  cracked  bodies.  Specifically,  a  laser  interferometry 
system  (including  optical  and  vibration  isolation  components  with  image  processing  hardware 
and  software),  and  a  tensile  substage  for  an  electron  microscope,  form  the  core  of  an  experimental 
program  which  complements  mature  analytical  and  computational  investigations. 

Significant  progress  has  already  been  made  on  two  particular  research  projects  made 
possible  by  this  grant.  The  appendix  includes  copies  of  publications  generated  from  the  initial 
stages  of  experimentation  in  both  investigations.  The  first  of  these  is  the  development  of  a 
coupled  experimental  and  computational  methodology  for  solving  inverse  fracture  problems.  This 
program  hinges  on  the  laser  holography  equipment.  The  methodology  can  be  utilized  m  the 
fracture  mechanics  of  homogeneous  and  heterogeneous  media  to  determine  the  location,  size,  and 
geometry  of  an  embedded  crack  given  the  measured  displacement  history  of  a  free  surface.  The 
technique  will  be  valuable  as  a  research  tool  for  those  interested  in  tracking  the  development  and 
propagation  of  cracks  inside  test  specimens.  The  second  program  involves  examining  microscale 
sliding  and  debonding  at  fiber-matrix  interfaces.  This  study  relies  upon  the  SEM  tensile  substage 
and  related  equipment. 

In  addition  to  these  initiatives  currently  underway,  the  equipment  is  flexible  enough  to 
provide  avenues  for  a  wide  variety  of  mechanics  studies.  For  example,  the  high  sampling  rate  of 
the  laser  holographic  system  will  permit  the  evaluation  of  the  dynamic  fracture  behavior  of 
composites,  the  nature  of  surface  cracks,  and  multi-site  damage  in  composite  panels.  The  1000  lb 
tensile  substage  for  the  SEM  will  allow  monitoring  of  the  micromechanical  behavior  of  a  wide 
range  of  materials,  including  metal  matrix  composites,  polymer  matrix  composites,  and  organic 
and  inorganic  structures  occurring  in  nature.  We  are  currently  completing  a  project  in  the 
biomimetics  of  fatigue  and  fracture. 


3 


Equipment  Acquired 


II.  Equipment  Acquired 

This  section  contains  a  list  of  the  research  equipment  purchased  with  the  funds  provided  by  grant 
number  F49620-93”  1-0550.  As  stated  in  the  letter  that  was  sent  to  you  on  July  1 1,1994  regarding 
an  extension  request,  the  Electronic  Speckle  Pattern  Interferometry  (ESPI)  system  that  was 
specified  in  the  proposal  was  no  longer  being  sold.  The  components  that  were  necessary  to  build  a 
similar  system  with  some  additional  features  are  presented  in  an  itemized  list  under  the  ESPI 
heading. 


ELH:T[RCNIC  speckle  pattern  INTERFERCMETRY  system  $41823 


TitM 

photographic  plate  holder 

holographic  plates 

rubber  tank  w/floating  cover 

optical  power  ureter  and  detector 

actuator  and  controller 

TI  123  CCD  camera 

25rnTi  lens 

computer 

L/3-3  frame  grabber  w/4  source  cable 
LabVIEW  software 
GPIB  beard 

ConceptVi  Level  2  imaging  software 

breadboard  (l'x2') 

support  posts  for  breadboard 

twD-stage  translation  stage 

translation  stage 

darrped  rod  w/rack 

rod  w/rack 

rod  clairp 

portable  air  compressor 

pressure  transducer 

digital  pressure  gauge 

top  beam  steering  accessory 

bottom  beam  steering  accessory 

variable  b^msplitter 

beamsplitter 

corpact  spatial  filter 

pinhole  aperture  (5  microns) 

elliptical  mirror 

Pyrex  mirror  (1"  dia) 

mirror  (2"  dia) 

mirror  nount 

minror/beamsplitter  mount 
objective  irount 
objective  (60x) 
riser  block 
base  plate 
rotatable  platform 
laboratory  jack 
screw  kit 


quanhiy 

3 

3 

3 

1 

1 

3 

3 
2 
1 
1 
1 
1 
1 
1 
1 
1 
2 
6 

4 
1 
1 
1 
4 
2 
1 
2 
2 
2 
7 
2 
2 
3 

3 
2 

4 
7 
7 
2 

5 
2 


VENDOR 

Jodon 

Jodon 

Calumet  Photographic 

Nayport 

Newport 

Symx) 

Syrreo 

Gateway  2000 
Scion 

National  Instruments 

National  Instruments 

G  T  F  S 

Newport 

Newport 

Edmund  Scientific 

Newport 

Newport 

Newport 

Newport 

Newport 

DCI  Instruments 

DCI  Instruments 

Newport 

Newport 

Newport 

Newport 

Newport 

Newpxart 

Newport 

Newport 

Newport 

Newport 

Newport 

Neport 

Newport 

Newport 

Newport 

Newport 

Newport 

Newport 
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Equipment  Acquired 


S^ILIZED  H^e  LASER  w/POWER  SUPPLY  $6584 


item 

VENDOR 

35irW  laser 

Spectra-Plysics 

VIBRATIOl  ISOLATICN  TABLE  w/SUPPORTS  $7954 

item 

VENDOR 

vibration  isolation  table 

Newport 

vibration  isolators 

Newport 

OPTICS  $1775 

item 

VENDOR 

zoom  microscope  w/base  and  post 

Edmund  Scientific 

Sory  C350  camera 

Edmund  Scientific 

rack  and  pinion  movement 

Edmund  Scientific 

avC  to  BNC  cable 

Edmund  Scientific 

1000  Ld  tensile  substage  $6300 


ITEM  VENDOR 

tensile  substage  w/extended  travel  Edmund  Scientific 


LOAD  READOUT  SYSTEM  $3095 


item  vemxr 

load  readout  systam  Ernest  Fullam 


MOTOR  DRIVE  SYSTEM  $2000 

item  vendor 

motor  drive  systan  Ernest  Fullam 
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AIAA-95-1507-CP 

Inverse  analysis  for  embedded  fractures 

William  D.  Keat*  and  Mark  J.  Hallidayt 
Clarkson  University 
Potsdam,  New  York 

Michael  C.  Larson^  and  Melody  L.  Arthur^ 

Tulane  University 
New  Orleans,  Louisiana 


Abstract 

A  coupled  computational  and  experimental 
investigation  was  undertaken  to  assess  the  feasibility  of 
a  procedure  for  subsurface  crack  identification  based  on 
inspection  and/or  inversion  of  surface  displacements. 
The  study  began  with  the  linear  problem  of  generating 
contour  maps  of  the  surface  deformations  produced  by 
buried  fractures  of  known  geometry  and  loading.  An 
indirect  boundary  element  formulation  using  the 
fundamental  solution  for  tensile  and  shear  multipoles 
near  a  half-space  provided  an  efficient  mathematical 
representation  of  the  3-D  fracture.  These  preliminary 
results  offered  evidence  for  the  existence  of  unique 
correspondences  between  crack  geometry  (and  loading) 
and  the  resulting  uplift  at  the  free  surface.  The  inverse 
problem  of  crack  identification  was  then  addressed 
beginning  with  the  development  of  a  hybrid  of  the 
Marquardt-Levenberg  algorithm.  Numerical  and  physical 
experiment  were  conducted  to  assess  robustness  of  the 
proposed  inversion  methodology.  The  experimental 
medium  was  a  cube  of  transparent  brittle  material  in 
which  a  fracture  was  hydraulically  pressurized. 
Displacements  induced  at  the  surface  of  the  specimen 
were  measured  by  laser  interferometry  and  compared  to 
numerical  results. 


Jntroduction 

This  paper  examines  the  feasibility  of  a 
nondestructive  evaluation  tool  based  on  using 
measurements  of  surface  deformations  to  estimate  the 
geomeuy  of  a  three-dimensional  embedded  fracture. 
Potential  applications  range  from  the  tracking  of 
underground  hydraulic  fractures  to  laboratory  monitoring 
of  crack  growth  in  test  specimens  and  the  nondestructive 

inspection  of  structural  components  under  in-service 
loadmg. 
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The  proposed  methodology  fits  the  standard 
definition  of  an  inverse  problem^  Uniqueness  of  the 
solution  is  therefore  a  central  issue.  In  fact,  there  are  an 
infinite  number  of  fracture  geometries  that  may  be 
associated  with  a  given  surface  deformation  pattern. 
Therefore  constraints  have  to  be  placed  on  the  set  of 
allowable  geometries  in  a  manner  consistent  with  the 
underlying  physics.  Furthermore,  any  computational 
strategy  must  be  robust  enough  to  deal  with  the  noise 
levels  arising  from  the  following  sources;  model 
approximations  of  material  behavior,  uncertainty  in  the 
boundary  conditions,  and  errors  in  the  experimental 
measurements.  The  other  issue  concerns  the  efficiency 
of  the  inversion  algorithm.  The  latter  must  consist  of 
two  parts:  a  model  of  the  three-dimensional  fracture  near 
a  free  surface,  and  a  nonlinear  equation  solver.  If  either 
one  proves  to  be  excessively  inefficient,  it  could  render 
the  completed  algorithm  impractical. 

Because  the  very  nature  of  the  inverse  problem 
relies  on  processing  measured  data,  surface 
displacements  in  this  case,  a  comprehensive 
methodology  must  synthesize  an  analysis  capability  for 
interpreting  measurements  with  an  experimental  data 
gathering  technique.  Practical  inputs  to  the  inverse 
solution  procedure  will  be  subject  to  the  stoichiometric 
considerations  which  necessarily  accompany  any 
metrological  endeavor.  Using  inputs  generated  by  a 
computer  code  to  validate  a  computational 
methodology’s  utility  as  a  practical  analysis  tool  would 
be  a  dubious  undertaking.  In  order  to  reveal  the  worth  of 
an  inverse  solver  as  a  useful  tool  for  researchers  and 
technicians,  one  must  demonstrate  the  accuracy  by 
inputting  actual  measurements  and  comparing  the 
results  with  known  solutions. 

The  following  paper  describes  progress  made 
towards  developing  a  coupled  - computational  and 
experimental  technique  for  solving  inverse  problems 
related  to  fracture  mechanics  of  homogeneous  and 
heterogeneous  materials. 


Fracture  Model 


To  be  suitable  for  the  inverse  problem,  a  fracture 
model  must  be  both  efficient  and  general.  For  example. 
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Mellings  and  Aliabadi^  employed  a  dual  boundary 
element  method  to  study  the  inverse  problem  in  two 
dimensions.  Those  endeavoring  to  solve  the  inverse 
problem  in  three  dimensions  have  mainly  focused  on 
deriving  appropriate  analytical  solutions.  By  necessity, 
fracture  models  of  this  type  tend  to  be  restricted  to  only 
the  simplest  of  geometries,  which  in  many  cases 
precludes  the  possibility  of  finding  a  solution  with  a 
high  level  of  confidence.  A  strategy  holding  far  greater 
promise  is  to  derive  the  displacement  influence 
functions  for  infinitesimal  fractnre  events  such  as  the 
multipole.  These  are  applicable  to  a  broader  range  of 
geometries  when  incorporated  in  an  integral  equation 
fonnnlation.  A  fairly  complete  summary  of  the  various 
3-D  solutions  can  be  found  in  a  paper  by  Okada^.  In 
some  instances  they  have  been  applied  to  investigate  the 
inverse  problem'*. 

The  fundamental  solution  for  point  forces  near  a 
planar  bimaterial  interface,  derived  by  Rongved^  in 
terms  of  Papkovitch  functions,  served  as  the  starting 
point  for  deriving  the  displacement  fundamental 
solution  for  tensile  and  shear  multipoles  in  a  semi¬ 
infinite  region^.  The  PapkovitcV  functions  were 
differentiated  twice,  first  with  respect  td  the  field  point 
coordinates  in  accordance  with  the  following  classical 
elasticity  equation  relating  the  Papkovitch  functions, 

B^'and  to  the  displacement  components,  : 


opening  and  shear 


^^3  )^i,j  | 


where  n,  J,  and  t  form  an  orthogonal  triad  of  local 
basis  vectors  oriented  with  n  normal  to  the  crack  plane; 
hP]  is  the  intensity  of  the  opening  mode  (or  tensile) 
multipole;  hP2  and  hPs  are  the  intensities  of  two  shear 
multipoles  which  are  superimposed  onto  the  tensile 
multipole  to  fully  represent  a  fracture  event  of 
infinitesimal  extent. 

The  displacements  induced  at  a  field  point  on  the 
free  surface  can  be  found  by  summing  the  effects  of  the 
multipoles  acting  over  the  crack  surface,  Sc-  If 
expressed  in  terms  of  components  of  relative 
displacement  5j  (i.e.  crack  opening,  dip,  Md  slip)  which 
are  related  to  the  multipole  intensities  through  Hooke's 
law,  the  following  integral  equation  results: 

(3) 


where  the  superscript  k  indicates  the  direction  of  the 
applied  force  and  v  is  Poisson’s  ratio.  A  second  set  of 
derivatives  was  taken  with  respect  to  the  source  point 
coordinates  consistent  with  multipole  definitions  for 


where  T?  are  the  displacement-influence  coefficients 
derived  in  equation  (2).  The  integral  is  nonsingular  for 
embedded  fractures  and  thus  can  be  evaluated  by  Gauss 
quadrature.  Note  that  the  key  feature  of  this  surface 
integral  method  is  that  only  the  surface  of  the  fracture 
needs  to  be  discretized. 

Demonstration  results  have  been  presented  in 
Figure  1  for  planar  and  nonplanar  fracture  geometries. 


Two  types  of  graphical  representations  have  been 
employed  to  represent  vertical  displacements  ( w)  at  the 
free  surface;  the  vertical  heights  of  points  lying  on  the 
'surface  plots’  are  proportional  to  the  value  of  w; 
'contour  plots'  present  the  same  results  as  lines  of 
constant  w .  Results  are  characterized  by  a  trough  (or 
depression)  along  a  line  approximating  where  the  kinked 
fracture  plane  intersects  the  free  surface.  In  general,  the 
trough  will  lie  directly  above  the  edge  of  the  aack  when 
the  fracture  is  perpendicular  to  the  free  surface,  and  will 
disappear  below  the  free  surface  as  the  parallel 
orientation  is  approached. 

Inversion  Algorithm 

The  determination  of  fracture  geometry  and 
location  from  experimental  measurements  of  surface 
displacement  is  formulated  as  a  nonlinear  least  squares 
problem.  This  is  a  logical  choice  since  inherent  errors, 
associated  with  the  experimental  observations  and  with 
restrictions  imposed  on  the  fracture  geometry,  serve  to 
make  an  exact  solution  highly  improbable.  The  scalar 
function  to  be  minimized  is  given  by: 

m  2 

R(x)  =  We^(y,.)]  (4) 

1=1 

where  x  is  the  vector  of  unknown  parameters  defining 
the  fracture  geometry;  {/?]  is  the  residual  vector;  w  and 
are  theoretical  and  experimental  values, 
respectively,  of  the  surface  displacement  component 
which  is  peipendicular  to  the  free  surface;  y,.  are 
coordinates  defining  the  location  of  the  /th 
measurement  site  on  the  free  surface.  Only  the  w 
component  of  surface  displacement  has  been  retained 
here  due  to  limitations  on  the  current  experimental 
setup. 

The  nature  of  the  proposed  application  required 
the  development  of  a  numerical  algorithm  with  the 
capability  to  solve  both  small  and  large  residual 
problems.  Following  in  the  footsteps  of  Dennis  et 
al.  ’  ,  a  classical  Marquardt-Levenberg  strategy  has  been 
combined  with  a  secant  approximation  to  secure  at  least 
superlinear  convergence  in  the  neighborhood  of  the 
global  minimum. 

To  begin  the  derivation  of  the  governing  matrix 
equation,  the  usual  statement  of  the  Newton  method  is 
rewntten  to  take  advantage  of  the  special  structure  of  the 
least-squares  equations: 

[jkJk  +  S^)Pk  =  -  JI  Rk  (5) 

where  Jf.  is  the  Jacobian  matrix,  Sf.  is  the  Hessian 
matrix,  p/.  is  the  search  vector  being  solved  for,  and  R^. 
denotes  the  residual  vector.  All  quantities  have  been 


defined  for  the  ^th  iteration.  At  the  conclusion  of  each 
iteration,  the  vector  of  unknown  fracture  parameters 
( )  is  updated  by  adding  in  the  search  vector  as  below: 

^k-¥\  =^k  +  Pk  (6) 

The  analytical  evaluation  of  the  derivatives 
required  for  and  is  all  but  impossible  due  to  the 

complexity  of  the  displacement-influence  functions  T? 

and  the  use  of  numerical  integration  procedures  with 
equation  (3).  Furthermore,  the  use  of  standard  finite 
difference  formulae  is  an  extraordinarily  expensive 
proposition  in  terms  of  CPU,  and  in  the  case  of  is 
prohibitively  expensive.  Thus  for  small  residual 
problems  it  is  both  'reasonable  and  advantageous  to 
neglect  S^. 

The  resulting  Gauss-Newton  formulation 
guarantees  a  descent  direction  but  does  not  assure  a 
function  decrease  as  over-shooting  may  occur.  To 
overcome  such  problems,  Levenberg?  proposed  the 
following  modified  form:  . 

{Jk  -^k  +  Pk^)  Pk  =  -  Jk  Rk  (7) 

where  I  is  the  identity  matrix  and  is  a  scaling 
factor.  As  the  value  of  Pi.  is  increased,  the  length  of  the 
search  vector  decreases  and  it  direction  approaches  the 
negative  gradient  (as  per  the  method  of  steepest  descent). 
Numerical  implementations  of  equation  (7)  have  been 
proposed;  the  most  popular  of  these  is  the  Marquardt- 
Levenberg  algorithm The  solution  of  equation  (7) 
was  carried  out  using  a  QR  algorithm  in  order  to  avoid 
the  numerical  instabilities  associated  with  the  normal 
equations  solution  of  small  residual  problems.  The 
least-squares  problem  is  therefore  cast  in  the  followins 
matrix  form: 


When  applied  to  small  residual  versions  of  the  fracture 
problem,  the  Marquardt-Levenberg  algorithm  exhibited 
the  expected  q-auadratic  convergence  and  was 
consistently  robust^. 

However,  with  large  residual  problems  the 
solution  progressed  at  such  a  slow  rate,  that 
modifications  seemed  in  order.  Mord’s^^  trust  region 
strategy  did  not  seem  to  provide  a  solution  because, 
though  robust,  it  is  still  essentially  a  small  residual 
algorithm.  The  best  of  the  remaining  choices  was  to 
adopt  one  of  several  available  secant  approximations  to 
The  original  Newton  formulation  is  once  again 
modified  to: 

{■^k-^k  +  A  +  Pti)  Pt  =  -  JI R^  (9) 
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where  Af.  represents  an  approximate  Sj^.  The 
recommendation  of  Dennis  et  al7>*,  regarding  the 
selection  of  A^,  is  to  choose  a  matrix  which  is 
sjTnmetric,  does  not  require  the  computation  of  any 
second  derivatives,  and  approaches  zero  as  the  residual 
approaches  zero. 

The  presence  of  A,^  in  equation  (9)  leads  to  some 
confusion  over  how  to  apply  the  QR  algorithm  in  order 
to  obtain  a  solution.  A  relationship  equivalent  to 
equation  (8)  can  be  established  but  it  requires  the 
Cholesky  decomposition  of  A,,.  This  does  not  appear  to 
be  a  viable  approach,  because  the  positive  definiteness 
of  Ai  has  not  been  enforced  and  a  perturbed 
decomposition  would  only  tend  to  corrupt  the  estimate 
of  St.  It  is  not  clear  how  Dennis  et  al.  (1981)  resolved 
this  issue.  The  solution  adopted  here  was  to  solve  the 
normal  equations  directly  using  LU  decomposition.  In 
the  event  that  the  residual  becomes  very  small, 
whereupon  the  normal  equations  solution  runs  the  risk 
of  going  unstable,  the  algorithm  switches  back  to 
equation  (7)  and  its  more  reliable  QR  solution 
procedure. 


Numerical  studies  were  conducted  to  assess  the 
range  of  applicability  of  the  inverse  algorithm. 
Existence  of  a  unique  global  minimum  was  promoted  in 
these  studies  by  constraining  tlie  fracture  to  be  planar, 
ellipUcal  in  shape,  and  of  constant  crack  opening.  In 
Keat  it  was  demonstrated  that  the  inversion  code  will 
converge  rapidly  for  both  large  and  small  residual 


problems  while  allowing  for  up  to  seven  degrees  of 
freedom  (e.g.  angles  of  orientation,  spatial  position, 
aspect  ratio  of  the  fracture,  crack  opening).  Recent 
emphasis  is  on:  (1)  defining,  as  precisely  as  possible, 
the  range  of  geometries  that  will  lead  to  convergence  to 
the  correct  global  minimum  and  (2)  the  sensitivity  of 
the  converged  solution  to  errors  in  both  the 
measurement  system  and  the  model  approximation. 

For  the  constrained  set  of  fracture  geometries 
under  consideration,  the  following  independent 
parameters  were  identified  as  having  a  significant  impact 
on  whether  or  not  convergence  would  occur:  number  of 
degrees  of  freedom,  nearness  of  the  starting  guess  to  the 
converged  solution,  nondimensional  crack  depth,  ratio 
of  crack  depth  to  laser-scanned  surface  width,  number  of 
measurement  sites  on  scanned  free  surface,  and 
systematic  and  random  errors  associated  with 
experimental  measurements. 

Sample  results  from  a  parametric  study  which 
was  undertaken  to  isolate  the  effects  of  each  of  these 
variables  are  given  in  Figures  2  and  3.  In  generating 
both  sets  of  results,  the  only  fracture-defining  variable 
that  has  been  allowed  to  change  is  the  crack  depth.  The 
perpendicular  orientation  of  the  test  geometry  with 
respect  to  the  planar  free  surface,  the  other  spatial 
coordinates,  the  radius  of  the  circular  crack  front,  and  the 
magnitude  of  crack  opening  have  aU  been  held  fixed. 

Figure  2  shows  how  bad  the  starting  guess  for 
depth  can  be  and  still  lead  to  convergence.  The  abscissa 
displays  the  crack  depth  corresponding  to  the  exact 
solution;  the  ordinate  gives  the  lower  limit  on  the  set  of 
starting  guesses  for  crack  depth  that  lead  to  conversence. 

It  is  interesting  to  note  that  the  ratio  of  the'’limit 
starting  guess  to  corresponding  actual  crack  depth 


10  15  20  25 

depth/radius 


Figure  2.  This  plot  shows  the  maximum 
depth  for  an  initial  trial  solution  that  will 
converge.  The  constant  variation  in  surface 
displacements  (bias)  is  given  in  percent 
of  the  maximum  dispacement  for  a 
reference  crack  with  a  nondimensional 
depth  of  2. 


10  15  ■  20  25  30 
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Figure  3.  This  plot  shows  the  error  in  the 
calculated  crack  depth  when  a  systematic 
bias  is  applied  to  the  surface  displacements. 
As  in  Figure  2,  the  constant  bias  is 
referenced  to  the  maximum  displacement 
for  a  crack  with  a  nondimensional  depth 
of  2. 
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remains  essentially  constant  over  the  range  of  depths 
examined.  Furthennore,  this  ratio  is  invariant  with 
respect  to  the  systematic  error  (or  bias)  for  small  biases. 
Not  shown  is  the  dependence  of  the  limit  starting  guess 
on  the  ratio  of  crack  depth  to  scanned  surface  wid^. 

Figure  3  shows  the  sensitivity  of  errors  in  the 
converged  solution  to  systematic  errors  associated  with 
the  field  point  measurements.  The  systematic  error  (or 
bias),  which  might  for  example  correspond  to  a  rigid 
body  movement  of  the  measured  object,  was  introduced 
by  adding  a  uniform  displacement  field  to  all  the  field 
pomt  measurements.  Percentage  biases  are  referenced  to 
the  largest  surface  displacement  produced  by  the  fracture 
model  at  a  nondimensional  crack  depth  of  2.  Results 
indicate  that  for  nondimensional  crack  depths  of  less 
than  20,  ^e  percentage  error  observed  in  the  converged 
solution  is  less  than  the  percentage  error  introduced  in 
the  surface  displacements. 


Laser  Interferometry  Measurement  .Svstpm 


The  preliminary  focus  of  the  experimental  phase  of  this 
project  was  to  induce  an  inlemal  fracture  of  known 
geometry  and  loading  within  a  medium  and  to  measure 
the  displacement  field  induced  on  one  of  its  free 
surfaces.  The  specimen  under  observation  was  a 
transparent  cube  of  polymethyl  metaciylate.  Cast  inside 
each  specimen  was  a  waxed  plastic  film  cut  to  the 
desired  crack  shape  which  forms  an  area  with  a 
predilection  for  fracture  formation. 

The  thin  aluminum  disc  was  mounted  on  the  end 
of  a  fluid  feed  tube  which  was  embedded  in  the  acrylic. 
A  penny-shaped  pre -crack  was  formed  by  attachina  a  6  4 
centimeter  diameter  piece  of  plastic  to  the  bottom'of  the 
diameter  aluminum  disc  using  tape.  In 
addition  to  providing  a  secure  hold  on  the  plastic  while 
me  acryhc  is  curing,  this  adhesive  tape  prevented  acrylic 
from  clogging  the  tube.  A  high  pressure  pump 
connected  to  the  fluid  feed  tube  supplied  intern^ 

Eert  IL  This  configuraUon 

allowed  the  predetermined  crack  to  be  unifoimly 

pressunzed  as  water  was  slowly  forced  down  the  feed 

The  out-of-plane  displacement  field  on  a  surface 

a  fracture  is  on  the  order  of 

magnitude  of  micrometers.  The  method  chosen  to 

h^ifa  \  out-of-plane  displacements  is 

^"ferferometry.  Whereas  devices  such  as 

avSace  dknlTr  extensometers  would  be  able  to  detect 
hnfn  °  ^/®P'2cements  over  finite,  discrete  lenaths 
to  interference  of  lig],; 

SiTce  ?uS  ""'ire 

resolve  "on-invasive  method  can 

of  nanometers  without 

ssoTrstsr"" 

emDloved^iL^‘?'''"'\P'"‘'  interferometric  technique 
P  yed,  two  laser  beams,  one  that  illuminates  the 


surface  of  the  specimen  and  one  that  serves  as  a 
reference,  interfere  with  one  another  on  a  recordinc 
medium  to  fonn  an  image  of  the  specimen.  When  twS 
such  images,  i.e.  one  of  a  specimen  in  its  undeformed 
state  and  one  in  its  deformed  state,  are  reconstructed  on 
the  same  hologram,  a  fringe  pattern  will  appear.  This 
patteni  is  indicative  of  the  amount  and  type  of  motion 
suffered  by  the  tes.  specimen  between"  e  fisTt  ^d 
second  exposure  of  the  holographic  medium. 

The  displacement  vector  of  an  arbitrary  point  on 
the  surface  of  the  object  is  a  simply  a  function  of  the 
optical  phase  difference,  0.  Considering  that  beams  of 
light  travel  in  wave  field  that  are  sinusoidal  in  nature 
the  optical  phase  difference  is  the  variance  in  phase  at 
which  two  beams  intersect  on  the  recording  medium  It 
follows  form  Huygeh's  Principle  of  Superposition  that 
when  two  light  beams  of  equal  magnitude  travelino  in 

f  ""k?  7^.  f  ^  interfere,  the  irradiancl  is 

oubled.  On  the  other  hand,  if  one  beam  lags  behind  the 
other  by  one  half  of  the  wavelength,  it  destructively 
mterferes  with  the  first  such  that  no  light  is  visible^ 
Quantitauvely,  the  phase  difference  between  the  beams' 
or  wave  fields,  is 


In 

T 


KL 


(10) 


where  X  is  the  wavelength  of  lioht  i  is  the 
displacement  vector,  K  is  the^ensitivity  vector,  and  k, 
and  are  unit  vectors  along  the  direcuon  of  viewing 
and  illumination,  respectively. 

The  optical  arrangement  for  evaluating  the 
surface  contours  of  the  acrylic  test  specimen  is  depicted 
in  Figure  4.  A  35  mW  helium  neon  laser  provided  light 
at  a  wavelength  of  632.8  nanometers.  The  li^ht  was 
divided  with  a  variable  beam  splitter  into  a  reference  and 
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Figure  4.  Holographic  apparatus  and  imase 
acquisition  system. 
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an  object  path.  Along  both  paths  the  beams  were 
reflected  and  spread  via  use  of  mirrors  and  spatial  filteis. 
These  beams  were  recombined  on  a  4  x  5  inch  glass 
plate  which  sat  in  a  microposidonable  plate  holder.  A 
capture  board  then  grabbed  the  images  viewed  through  a 
CCD  camera  and  transferred  them  to  a  computer*^ for 
image  processing. 

The  viewing  and  illumination  vectors  were 
perpendicular  to  one  another  and  were  oriented  such  that 
the  surface  normal  of  the  specimen  bisects  them.  Thus 
the  sensitivity  vector,  K,  was  reduced  to  single 
component  noimal  to  the  surface. 

The  surface  of  the  acrylic  specimen  was  painted 
white  to  reflect  diffusely  the  laser  light  through  the 
holographic  plate  and  create  an  original  image  of  the 
surface  that  appears  bright.  As  the  object  deformed,  the 
light  reflected  from  each  point  on  its  surface  traveled  a 
slightly  different  distance  than  when  taking  the  orisinal 
image.  If  the  distance  was  altered  by  half  the  wavelensth 
of  the  laser  light,  this  corresponding  point  viewTed 
through  he  hologram  will  change  from  bright  to  dark. 
Substitution  of  n  radians  (one  half  wavelength)  for  ^ 

In  equation  (10)  suggests  that  neighboring  dark  fringes 
have  a  difference  in  ou^of-plane  displacements  of  224 
nanometers. 

The  pre-cracked  specimen  and  the  optical 
components  were  secured  to  a  vibration  isolation  table 
to  eliminate  erroneous  fringes  stemming  from 
background  excitations.  A  specimen  holder  ensured  that 
there  would  be  no  rigid  body  motion  durino 
pressurization.  * 

Before  supplying  pressure  to  the  test  specimen,  a 
holographic  exposure  of  the  original  surface  was  taken 
and  the  plate  was  developed.  After  carefully 
reposmonmg  the  plate  in  its  holder,  the  CCD  camera 
began  transferring  images  at  thirty  frames  per  second  to 
the  computer  bus.  The  first  frame  captured  contained  tlie 


holographic  image  of  the  specimen  in  its  undeformed 
state;  the  ensuing  frames  were  holographic  images  of 
the  specimen  at  subsequent  levels  of  internal 
pressurization.  The  camera  was  oriented  such  that  each 
picture  contained  a  holographic  image  of  the  defonnina 
surface  as  well  as  a  digital  display  of  the  corresponding 
pressure  applied  to  the  pre-cracked  region.  The  procedure 
allowed  the  pressure  to  the  specimen  to  be  steadily 
maeased  while  recording  a  history  of  the  fringes  as  they 
develop  and  evolve  into  different  shapes.  Monitoring 
fringes  as  they  originate,  relocate,  or  disappear  from  the 
surface  provides  a  means  for  tracking  the  order  of  the 
fringes. 

To  evaluate  a  displacement  contour  of  the  surface 
for  a  chosen  internal  pressure,  the  displacement  field 
corresponding  to  the  initial,  unpressuiized  holographic 
image  of  the  test  specimen  was  subtracted.  Because  the 
plate  holder  is  only  micropositionable,  realigninc  the 
developed  plate  to  a  sensitivity  in  the  submicron  region 
is  not  achievable.  Therefore,  while  readjustments  to  the 
plate's  position  in  the  holder  was  made  to  eliminate 
^ost  all  of  the  fringes,  a  few  remained.  Fringes 
induced  by  pressuriz§tion  evolved  from  the  inifial 
pattern.  Tracking  fringes  from  the  image  corresponding 
to  the  undeformed  state  to  its  new  position  at  some 
applied  pressure  yielded  the  developing  deformation  field 
of  the  specimen  surface. 


Preliminary  Experimental  Findino'; 

Figure  5  shows  a  sample  comparison  between  an 
experimentally  determined  out-of-plane  displacement 
field  with  a  computationally  generated  prediction.  This 
particular  correlation  is  precursory  in  that  it  was 
conducted  to  ascertain  whether  the  computational 
algorithm  generates  results  which  are  compatible  with 


those  yielded  from  the  experimental  measurements. 

The  experimental  field  is  noticeably 
asymmetrical  for  this  test  case.  This  is  tied  to  the 
manner  in  which  the  pre-crack  was  embedded  into  the 
specimen.  As  indicated  earlier,  a  thin  waxed  film  was 
attached  to  the  disk  on  the  end  of  the  fluid  feed  tube 
using  double-sided  tape.  The  tape  did  not  release 
uniformly  during  pressurization  and  hence  led  to  an 
asymetrical  crack  opening  distribution  within  the 
fracture.  This  problem  is  being  remedied  by  eliminating 
the  film  and  forming  the  existing  disk  to  the  desired 
geometry. 

In  spite  of  this,  one  can  clearly  see  a  decided 
qualitative  agreement  between  the  patterns,  especially  in 
the  depression  centered  on  the  crack  plane  flanked  by 
steep  outward  crests.  This  comparison  is  particularly 
encouraging  in  light  of  the  approximate  nature  of  the 
computational  model  employed.  The  fundamental 
solutions  for  a  semi-infinite  half  space  were  used  which 
only  account  for  one  of  the  free  surfaces  of  the  cubic 
test  specimen.  The  quantitative  correspondence  is 
expected  to  improve  dramatically  through  linking  the 
surface  integral  representation  of  the  fracture  with  a 
finite  element  formulation.  This  will  permit  the  explicit 
representation  of  the  bounded  geometry. 


Conclusions 

The  results  outlined  here  for  stationary  cracks 
lend  o-edence  to  the  strong  possibility  for  developing  a 
practical  nonlinear  inversion  methodology  for  resolving 
subsurface  damage  from  nondestructive  surface 
measirements.  The  numerical  evaluations  of  solution 
sensitivity  are  continuing  and  are  guiding  refinements  in 
the  experimental  measurement  technique.  The 
experimental  program  is  expanding  to  include 
observations  of  evolving  displacement  patterns 
associated  with  propagating  fractures.  This  history 
dependence  is  being  used  to  remedy  some  of  the 
uniqueness  problems  associated  with  partially  defined 
continuum  problems. 
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ABSTRACT 

Increasing  the  toughness  of  brittle  ceramics  with  fibers  has 
received  much  attention  in  the  last  five  years.  Many  models 
have  been  put  forth  to  describe  the  various  facets  of  cracking 
in  these  materials.  The  intent  here  is  to  evaluate  how  well  a 
relatively  simple  computational  representation  of  the  crack 
openings  for  a  fracture  in  a  fibrous  material  predicts  the  crack 
openings  actually  measured  in  an  experimental  specimen.  The 
cracks  are  modeled  as  distributions  of  discontinuities  within  a 
surface  integral  representation.  The  influence  of  the  fibers  are 
accounted  for  through  closing  tractions  applied  at  the  fracture 
surfaces.  Experimental  specimens  were  loaded  using  a  tensile 
sub  stage  mounted  in  a  scanning  electron  microscope  to 
visually  monitor  the  crack  and  damage  development.  Some  of 
the  physical  processes  associated  with  matrix  crack-fiber 
interactions  were  also  evident. 


INTRODUCTION 

A  major  driving  force  for  developing  ceramic  matrix 
composites  is  the  potential  of  fabricating  materials  which 
retain  the  excellent  thermal  and  stiffness  characteristics  of 
monolithic  ceramics  but  have  enhanced  fracture  toughness 
(e.g.,  Chou  et  al.,  1986).  Such  materials  would  find  obvious 
application  in  high  temperature  engine  components.  There  are 
a  number  of  mechanisms  at  the  microstructural  level  which 
hold  promise  for  achieving  this  end,  one  of  them  being  fiber 
reinforcement. 

The  fracture  of  fiber-reinforced  ceramics  is  typically 
different  from  that  in  monolithic  ceramics,  which  often  fail  by 
the  growth  of  a  single  crack  on  a  plane  normal  to  the  maximum 
principle  stress.  Fiber  composites,  on  the  other  hand,  can  fail 
by  a  variety  of  mechanisms,  depending  on  the  fiber/matrix 
interfacial  characteristics  and  the  interaction  between  a  crack 
front  and  the  reinforcement.  Crack  deflection,  fiber  pull-out, 
matrix  microcracking,  crack  bowing  and  fiber  bridging  have 
been  identified  as  major  mechanisms  of  toughening  of  fiber- 
reinforced  ceramic  matrix  composites. 


This  work  establishes  a  two-dimensional  computational 
model  to  represent  the  increased  toughness  exhibited  by  a 
[(0/90)4]s  composite  system  composed  of  a  glass  matrix 
reinforced  with  nicalon  fibers,  subjected  to  uniaxial  tension 
parallel  to  the  0*  fiber  direction,  through  crack  opening 
displacement  calculations.  This  approach  is  relatively  simple 
and  flexible  to  implement  in  comparison  to  many  other  two- 
dimensional  models  which  have  been  proposed  in  the  literature 
(Evans,  1990).  A  number  of  approaches  have  been  proposed  in 
the  literature  to  explain  the  increased  toughness  exhibited  by 
fiber  reinforced  ceramic  matrix  composites.  These  approaches 
are  based  upon  various  descriptions  of  how  to  account  for  the 
presence  of  the  fibers.  Common  models  include  shear  lag  ,  load 
transfer,  debonding,  and  the  fiber  pull-out  model.  The 
enhanced  toughness  of  fiber-reinforced  ceramic  composites  is 
believed  primarily  to  stem  from  the  fibers  left  intact  in  the 
wake  of  an  advancing  fracture.  This  phenomenon  is  known  as 
'bridging.’  In  the  above  models  it  is  common  to  represent  the 
bridging  fibers  as  closing  tractions  appled  to  the  crack  faces 
(e.g..  Marshal,  et  al.,  1985,  Marshal  and  Cox,  1987).  The 
intent  here  is  to  see  if  such  a  simplified  approach  captures  the 
proper  crack  opening  behavior  in  a  model  composite. 

COMPUTATIONAL  APPROACH 

Computational  implementation  of  the  two-dimensional 
displacement  discontinuity  method  for  fractures  was 
popularized  by  Crouch  (1976).  The  methodology  involves 
placing  displacement  discontinuities  along  the  boundaries  of 
the  region  to  be  analyzed  (in  this  case  a  crack),  then  solving  a 
system  of  algebraic  equations  to  find  the  discontinuity  values 
that  produce  prescribed  boundary  tractions.  The  problem  of 
determining  the  displacement  fields  in  an  elastic  half-space 
containing  an  edge  crack  is  of  particular  significance  in  this 
investigation  since  it  matches  the  experiments  conducted  for 
verification.  The  following  section  briefly  describes  the  two- 
dimensional  (plane  strain)  formulation  of  the  displacement 
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discontinuity  method  for  homogeneous,  isotropic,  linearly 
elastic  materials  for  semi-infinite  regions. 

The  analytical  formulation  for  the  half  space  is  obtained  by 

J  method  of  images.  Using  this 

meth^,  the  solution  is  found  in  two  stages.  In  the  first  stage 
of  tte  analysis,  an  mfinite  body  containing  two  fine  segments 

Sir  Hit  ’  m  y  >  0-  represents  the 

“.wt-  '^h.le  the  other,  in  y  <  0,  represents  its 

image  ,  reflected  about  the  line  y=0.  This  makes  the  shear 
action  _  <7^  yanish  all  along  y  =  0  due  to  symmetry.  It 

toe  t  -  n”  tZ  ^'’se  crack  Ues  along  ihe 

tn  t  ^  supplementary  solution  is 

added  to  make  the  normal  traction  /^  =  <7^  on  y  =  0  yanish. 

/loS  Z  deyelopment  as  Resented  in  Crouch 

(1^76)  with  appropnate  specialization.  Coordinates  |,  C  are 
related  to  the  points  in  the  x-y  system  by 

^  =  ri(x-b)-Hy-c) 
f  =  l(x-b)+n(y-c) 

m  which  I  and  n  are  direction  cosines  of  the  C-axis 
relatiye  to  the  x-  and  y-axes  and  b  and  c  are  the  x  and  y 
distances  to  the  origin  of  the  system.  The  image 

coordinates  ^  ^  are  described  by  the  transformation 


I  =  n(x-b)+l(y+c) 

C  -  ~Kx-b)+n(y  +  c) 


(2) 


^Denoting  the  stresses  due  to  the  actual  displacement  by 
0,j,  those  due  to  its  image  by  cr'  and  those  resulting  from  the 

supplementary  solution  by  of ,  the  complete  solution  for  the 
half -plane  y  >  0  is: 


<^:-=  cf^+al+a?, 
‘J  y  y  II 


(3) 


**  modeled  as  a  series  of  these  discontinuities 
distnbuted  oyer  discrete  elements.  The  discontinuity  is  diyided 

element  on  the  final  result  will  be  discussed  later  Also  the 
md^mts  of  these  elements  are  taken  as  the  points  wtorl  Se 

to  hf  "!  T  displacements  S  need 

Sn  h  ^  f  ^  m‘®"sity  factor 

M  be  calculated.  In  the  present  analysis,  S  is  taken  to  be  a 

rrem7nfr  7""  tT' 

1’  2  ’  ^3  for  each  element,  we  seek  to 

ormulate  the  problem  as  a  system  of  algebraic  equations: 

(4) 

In  this  equation,  C^are  influence  coefficients  calculated  using 
an^ytical  stress  relationships  and  r,  are  the  tractions  at  the 

.e«n.  load  bei.* 

‘‘  (5) 


where  n.  is  the  normal  to  each  element,  defining  the  traction 
direcfion.  The  most  computationally  intensiye  phase  of  the 
sol^ution  IS  calculating  the  nonsymmetric,  fully  populated 
influence  coefficient  matrix.  ^  ^ 

■  “umerical  procedure  inyolves  considering  the  semi- 
infinite  body  containing  N  crack  segments  across  which  the 
normal  displacements  are  discontinuous.  Let  the  jtb  segment 
have  half-width  a.,  direction  cosines  Ij  andrtj,  and  midpoint 
crordinates  Xj  ,  y.  and  the  normal  displacement  discontinuity 
Sj .  The  stress  <7,^  at  an  arbitrary  point  x,  y  in  the  body  due  to 
toe  ^  discontinuity  can  be  computed  directly  from  equation 
tp'™  *  total  stress  at  point  x,  y  is  given  by  a  summation  of  N 
ims,  one  for  each -.of  the  elemental  discontinuities  in  the 

c^a  ^  normal  traction  at  toe  midpoint  of  the  ith 

segment  can  be  put  m  the  form 


Kl- 


N 

s 

j  =  l 


(Sm)y  '*'(Sm),yj 


(6) 


;nr^  coefficients  in  the  above  expression  are  toe 

infinite  body  fiincUons  whUe  the  primed  coefficients  represent 
I^IIe  account  for  the  presence  of  toe  traction-free 

frp/ll  f**  present  case  of  an  edge  crack  perpendicular  to  toe 
free  surface,  aU  the  crack  segments  are  placed  end  to  end  and 

sIItem'^F!!  with  respect  to  the  global  coordinate 

system,  ror  this  case, 


(%), "  - 


A 


1 


1 


1 


2;r(l-v)  (^  +  a)  (4-a)  (^'  +  a) 
—5 - +{3y+3(^-c))[+ - 1_. 

1 


(7) 


rJ+2y(^-c)[' 


-] 


where  ^  and  are  given  by  equations  8  and  9. 

5,-  =  nj(X.-X.)-l.{y^-y.) 

=  h i^i-Xj)+ Tij (y, -yj) . 

=  nj{x,-Xj)+lj{y.+y.) 

Equation  (7)  is  toe  final  form  of  the  equation  used  to 
alculate  the  influence  coefficients  in  the  program  module, 
or  linear  elastic  bodies,  toe  stress  intensity  factor  K  is 

propensity  for  crack  extension:  the  higher  is 
the  stress  intensity  factor,  the  greater  is  the  tendency  for  the 
crack  to  advance.  Although  it  is  understood  that  toe  presence 
of  fibers  will  introduce  nonlinear  behavior,  toe  approach  taken 


(8) 


(9) 
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here  for  assessing  the  influence  of  the  bridging  fibers  is  to  use 
the  crack  opening  displacements  obtained  from  the  procedure 
outlined  above  to  calculate  an  effective  stress  intensity  factor. 
For  this  we  employ  the  asymptotic,  near-tip  crack  opening 
displacement  solution 


length  of  the  crack  behind  the  bridging  zone  over  which  the 
fibers  break. 

According  to  Budiansky  and  Amazigo  (1989),  the  crack 
opening  displacement  in  terms  of  the  axial  stress  in  the  fiber 
is  given  by 


S(r)  = 


K^^rK 


(10) 


where  k  =  3— 4v  for  plane  strain  and  k  =  -  for 

1  + V 

plane  stress,  G  is  the  shear  modulus,  r  is  the  distance  from  the 
crack  tip  to  the  center  of  the  crack  element  and  6(r)  is  the 
crack  opening  displacement. 

In  equation  (8),  r,  and  the  corresponding  S{r)  need  to  be 
chosen  as  close  to  the  crack  tip  as  practical.  It  is  possible  to 
expliciUy  incorporate  the  form  of  this  asymptotic  solution 
into  the  interpolation  function  of  the  elemental  displacement 
discontinuities  to  form  specialized  crack  tip  elements.  While 
this  is  advisable  for  monolithic  bodies,  we  do  not  take  this 
approach  here  since  we  do  not  know  how  the  non-linearities 
associated  with  the  fibers  and  interfaces  in  the  crack  tip 
process  zone  will  affect  the  nature  of  the  opening 
displacements.  Hence,  collocation  point  displacements  are 
fitted  to  match  the  parabolic  opening  near  the  crack  tip. 

The  computational  approach  was  checked  against  the 
theoretical  solution  for  a  center  crack  in  an  infinite  plate  to 
determine  heuristic  rules  for  extracting  the  stress  intensity 
factor  using  constant  elements.  A  trade-off  must  be  made 
given  that  (i)  the  constant  elements  do  not  capture  the  large 
deformation  gradient  near  the  tip,  so  the  more  of  the  crack 
which  is  considered,  the  better  the  actual  opening  is 
represented,  and  (ii)  the  asymptotic  crack  opening  solution 
used  to  calculate  the  stress  intensity  factor  breaks  down  as 
distance  from  the  tip  increases.  Figure  3  shows  a  plot  between 
the  calculated  stress  intensity  factor,  Kcal,  and  the  theoretical 
value,  Kth>  versus  percentage  of  the  crack  used  in  the  Kcal 
approximation.  Solutions  run  for  increasing  numbers  of 
elements,  i.e.  decreasing  element  sizes,  suggest  that  the 
portion  of  the  crack  to  include  in  the  stress  intensity  factor 
calculation  is  between  10%  and  20%  of  the  half  crack  length. 

Concerning  the  closing  traction  exerted  by  individual 
bridging  fibers,  it  is  likely  that  in  actuality  the  magnitude 
increases  with  increasing  opening  displacement  until  the 
fibers  break.  From  that  point  onwards  the  closing  traction 
decreases  until  the  fibers  are  completely  pulled-out  of  the 
matrix.  In  the  present  model  proposed,  this  cycle  is 
represented  by  a  constant  closing  traction  as  an 
approximation.  This  closing  traction,  T,  exerted  by  the  fibers 
is  applied  across  the  fiber  sites  (Figure  4)  and  is  superimposed 
over  the  opening  traction,  P.  The  relative  size  of  the  elements 
which  have  the  closing  traction  applied  is  tied  to  the  volume 
fraction  of  fibers  and  to  the  shear  resistance  of  the  interface. 

The  effect  of  the  bridging  length  and  the  closing  tractions 
exerted  by  the  fibers  on  the  effective  stress  intensity  factor  of 
^e  composite  is  studied.  The  computational  analysis  has  been 
extended  to  encompass  varying  volume  fractions  of  the  fibers 
in  the  composite.  In  the  present  model,  a  steady  state 
assumption  has  been  made  regarding  the  bridging  length.  This 
means  that  the  length  by  which  the  crack  extends  equals  the 


(11) 


where  is  the  axial  stress  in  the  fiber,  <7;.  is  the  reference 
stress,/  is  the  volume  fraction  of  the  fibers,  and  is  the 
critical  matrix  energy-release  rate. 

The  reference  stress  is  given  by 


1/4 

1 _ _ 

(12) 


where  is  Poisson's  ratio  for  the  composite,  E  is  the 
composite  modulus  of  elasticity,  and  E^  are  the  Young’s 

Moduli  of  the  fiber  and  matrix  respectively,  and  r  is  the 
radius  of  the  fiber. 

B  is  a  constant  and  is  given  by 


B  = 


2  (1-/)^ 


|l/4 


[6  log/+3(l>/)(3-/) 


(13) 


Equation  (11)  can  be  used  to  calculate  5'  by  substituting 
<J^  for  <7^,  where  3'  is  crack  opening  displacement  at  the 
back  of  the  bridge  zone  and  cr^  is  the  breaking  stress  of  the 
fiber.  This  gives  an  estimate  of  the  bridge  length,  when 
comparing  3'  with  the  crack  opening  displacements,  over 
which  the  closing  tractions  will  be  applied. 

EXPERIMENTAL  DETAILS 

The  substage  used  in  the  experiments  permits  a  user  to 
observe  the  effects  of  tensile  or  compressive  stresses  of  up  to 
1(X)0  lb  on  a  relatively  small  test  specimen  mounted  in  an 
SEM  chamber.  This  substage  has  a  load  cell  and  a  motor  drive 
system. 

T^e  specimens  are  made  of  a  composite  system  which 
consists  of  a  barium-magnesium  aluminosilicate  (BMAS) 
matrix  that  contains  SiC  over  BN  coated  Nicalon  fibers.  The 
fiber  coatings  have  been  applied  by  a  chemical  vapor 
deposition  process.  Carbon  is  deliberately  added  to  the  BN 
precursor  since  it  has  been  found  that  the  carbon  addition 
minimizes  reactivity  during  deposition  between  the  BN  and  the 
Nicalon  fibers. 

Experimental  specimens  were  cut  from  a  10  x  10  cm 
[(0/90)4]s  BMAS  glass-ceramic  composite  panel.  The  material 
was  hot-pressed  subject  to  a  maximum  temperature  of  1420^  - 
1440  C  at  6.9  MPa  pressure  for  5  to  10  minutes  in  an  Argon 
atmosphere.  Complying  with  the  requirements  of  the  tensile 
substage,  the  dimensions  of  the  experimental  specimens  was 
chosen  to  be  3.8  cm  long  x  0.65  wide  x  0.32  cm  thick.  A  small 
notch  was  machined  into  each  specimen  in  order  to  localize  the 
damage  and  crack  development. 

During  loading,  all  five  of  the  specimens  tested  exhibited 
the  same  stages  of  damage  evolution.  Figure  5  schematically 
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T^LE  1.  PREDICTIONS  AND  MEASUREMENTS 
FOR  CRACK  OPENING  DISPLACEMENTS  (COD) 


Distance 
Behind  The 
Crack  Tip  (cm) 


COD 

(numerical) 


COD 

(experimental) 


0.001 

1.1x10-6 

V  7 

0.91x10-6 

0.005 

2.0x10-6 

1.80x10-6 

0.010 

3.1x10-6 

2.70x10-6 

0.017 

3.8x10-6 

3.60x10-6 

0.020 

4.1x10-6 

4.50x10-6 

shows  the  development  of  failure  observed  in  the  notch  region 

involves  the  formation  of 
distnbuted  damage  in  the  matrix  around  the  notch.  Stage  two 
IS  characterized  by  the  formation  of  (usually  two)  microCTacks 
emanating  the  notch.  Finally,  in  stage  three,  one  of  the 
microscopic  cracks  grows  into  a  macroscopic  one.  This 
mcrocrack  extends  in  a  quasi-stable  manner  for  a  period  of 
time  before  dynamically  traversing  the  specimen.  During  the 
quasi-stable  extension  the  mechanisms  of  crack  bridging  and 
fiber  puU-out  can  clearly  be  seen  under  the  SEM.  Figme  6 
shows  a  photograph  of  the  trailing  edge  of  a  bridging  zone 
capturmg  the  transition  from  bridging  to  pull-out. 

COMPARISON  WITH  EXPERIMENT 

lie  results  of  the  numerical  model  are  compared  with  the 
expenmental  results  for  one  of  the  specimens  tested.  During 
e  quasi-stable  extension  of  the  dominant  micro-crack  the 
crack  openmg  displacements  were  recorded  under  the  scanning 
electron  microscope  under  a  load  of  880  N.  At  the  instant 

c^sen,  the  crack  was  2.5x10-4  m  long  and  subject  to  an 
openmg  pressure  of  44  MPa.  j  ^ 

®  ^  concentration  around  the  notch  region  and 

hence  the  actual  stress  experienced  by  the  crack  is  greater  than 
a  applied  remotely  to  the  specimen,  CTo.  The  stress 
concentration  effect  of  the  notch  is  approximated  by 

'  ‘’•(‘"f]  (14) 

where  is  the  stress  in  the  notch  region,  a  is  the  semi- 
ellipse,  and  b  is  semi-minor  axis  of  the 

Hence,  is  the  opening  traction  used  in  the  calculations. 
Also,  a  nominal  closing  traction  equal  to  is  assumed  in  the 
numencal  model.  The  major  and  minor  axes  of  the  elliptic 
notch  was  measured  and  equation  12  used  to  calculate  <t  c 
was  found  to  be  440  MPa,  ten  times  the  applied  nomind 

rrarV  H  through  13  are  used  to  calculate  S',  the 

crack  opening  displacement  at  the  back  of  the  bridged  zone. 

®i"S  -  15  GPa,  Ej  =  250  GPa,  =  96  Nm'l  ,  =2.5 

GPa,-/  =0.5,  r=5/mi,  E  =  132  GPa,  =  0.2,  5'  comes  out 

when  estimate  of  the  bridge  length 

unreinfor?E|e  ®  displacements  of  the 


rrSIV  ^^eck  opening  displacements 

(COD)  for  toe  edge  crack  in  a  half  plane  and  perpendicular  to 
the  free  surface  as  a  function  of  the  distance  behind  the  crack 

tip. 

Comparison  of  the  experimental  crack  opening 
displacements  with  the  numerical  crack  opening  displacements 

‘difference  ranges  from  about 
5%  to  20%.  It  can  be  seen  from  Table  1  that  the  numerical 
model  overestimates  the  crack  opening  displacements  in 
companson  to  the  experimental  crack  opening  displacements. 
The  major  problem  encountered  in  interpreting  the  measured 
data  was  toat  the  crack  trajectory,  although  straight  on  a 
macroscopic  scale,  appears  to  meander  across  the  specimen 
under  such  high  magnifications.  This  makes  it 


CONCLUSIONS 

The  metiiodology  can  be  used  to  predict  the  variations  in 
crack  openmg  profiles  with  the  volume  fraction  of  and  the 
frictional  resistance  afforded  by  the  fibers.  This  variation  is 
presumed  to  dictate  the  propensity  of  fracture  growth  and  can 
therefore  be  considered  an  "effective"  stress  intensity  factor, 
K^.  This  value  is  used  for  comparision  purposes  only  as 
being  indicative  of  the  severity  in  crack  opening  profiles 
Figure  7  shows  the  ratio  of  the  effective  stress  intensity 
factor  and  that  which  would  occur  in  an  unreinforced  material 
possessing  no  fibers  but  having  the  average  composite 
p  perries,  KcffiK^o  fiber »  versus  fiber  closing  traction,  T/P 
(refer  to  figure  4)  for  an  edge  crack  in  a  material  system  with 
vapung  volume  fractions  of  fibers.  Figure  8  shows  the  plot 
oetween  KeffKno  fiber  versus  //a  for  varying  T/P  ratios  for  an 
edge  crack,  where  //a  is  the  bridge  length  normalized  by  the 
crack  length.  Knowing  the  closing  traction,  we  can  find  the 
effective  stress  intensity  factor  from  Figure  7  for  a  given 
volume  fraction  of  fibers.  From  the  calculated  effective  stress 
intensity  factor,  the  bridging  length  can  be  estimated  from 
Figure  8  for  the  given  T/P  ratio. 

r  “  concluded  that  the  present  simple  approach  for 
mcjdehng  the  fiber  reinforced  ceramic-matrix  composites 
holds  promise  for  accurately  representing  crack  displacements. 
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FIGURE  1.  COORDINATE  SYSTEMS  FOR  ACTUAL  AND  IMAGE  DISPLACEMENT  DISCONTINUITIES, 


FIGURE  2.  AN  EDGE 


CRACK  SUBJECT  TO  MODE  I  TENSION.  THIS  GEOMETRY  MATCHES 


THE  EXPERIMENTS. 
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Iengtl3,^l^ack 

length 


LENGTH.  (N  IS  THE  TOTAL  NUMBER  OF  ELEMENTS) 
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distributed  damage  forms 
in  the  notch  region 


microcracks  emanate  from 
the  notch  at  about  45*  off-  axis 


one  microcrck  develops  into 
a  dominant  fracture 


FIGURE  5.  EVOLUTION  OF  DAMAGE  AND  CRACKING  OBSERVED  IN  THE  NOTCH  REGION 

OF  THE  EXPERIMENTAL  SPECIMENS. 


figure  6.  REPRESENTATIVE  SEM  PHOTOGRAPH 


OF  THE  END  OF  THE  BRIDGING  ZONE 
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FIGURE  7.  COMPUTATIONAL  PREDICTION  OF  THE  VARIATION  IN  NORMALIZED  STRESS  INTENSITY  FACTOR 
WITH  THE  FRACTION  OF  THE  CRACK  LENGTH  BRIDGED  BY  FIBERS  FOR  THREE  VALUE^F 
CLOSING  TRACTION.  ALL  THE  CURVES  REPRESENT  50%  VOLUME  FOACTraN 
OF  FIBERS  IN  THE  NICALON-GLASS  COMPOSITE. 


FIGURE  8.  COMPUTATIONAL  PREDICTION  OF  THE  VARIATION  IN  NORMALIZED  RTRFEc;  imtrmqitv 

^^tion  for  four  values  of  fiber  “LiMEF^rira^ 
REPRESENTSBR®G,NGLFNGT^^^^^^^^^^^ 

O  THE  FIBER  BREAKING  STRESS  FOR  THE  NICALON-GLASS  COMPOSITE. 
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